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Abstract 

We show that solitary waves for the 2D Euler-Korteweg model for capillary flu¬ 
ids display nonlinear instability when subjected to transverse perturbations. 
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1 Introduction 

We consider the motion of a compressible, inviscid and isentropic planar fluid, in which 
internal capillarity is taken into account. This phenomenon occurs for example at diffuse 
interfaces in liquid-vapour mixes |1]. In this model, the free energy of the fluid depends 
on both the density of the fluid, the scalar function p, and its gradient Vp in the following 
way: 

F(ftVp) = ri(p) + y(p)|Vp|^ 

with K and Fq two given smooth, positive functions. We then derive the pressure from 
the free energy as follows: 

r)F 1 

P(p, Vp) = p— - F = Po(p) + ^{pK\p) - V(p))| Vpp, 

in which Pq is the standard part of the pressure. Let Po(p) be the bulk chemical potential 
of the fluid, so that pp^p) = Po(p)- Then, the principles of classical mechanics yield the 
Euler-Korteweg equation that we will study: 

{ dtp + div {pu) = 0 

dtU + {u-V)u = V(V(p)Ap + lV'(p)|Vp|2-po(p)) (1) 

{p,u)\t=o = (po,Mo)- 

*Sorbonne Universites, UPMC Univ Paris 06, UMR 7598, Laboratoire Jacques-Louis Lions, F-75005, 
Paris, France 

'I'CNRS, UMR 7598, Laboratoire Jacques-Louis Lions, F-75005, Paris, France 


1 





The variables are t G M"*" and {x,y) G as is standard, the operators V, div and A 
contain only derivatives with respect to the space variables x and y. The unknowns of 
equation ([1]) are the density p and the velocity vector held u : M"*" x ^ R^. The scalar 
functions Qq and K are given. 


In this paper, we will be interested in the transverse stability of solitary wave solutions 
of dH). These are ID travelling waves written as 


Qc{t,x) 


pcit.x) 

Uc(t,x) 


Qc{x - ct), 


with Uc scalar (not a 2D vector held). Based on a remark by T. Benjamin [T], S. Benzoni, 
R. Danchin, S. Descombes and D. Jamet showed in [1] that the hamiltonian structure 
of the system led to the existence of solitary wave solutions for every c G R in the case 
of a van der Waals-type pressure law. There are two kinds of travelling-wave solutions, 
depending on the homoclinic or heteroclinic nature of the connecting orbits. In the 
homoclinic case, the wave has identical endstates, and we can write 


lim Qc('^) Qc ,00 (PcjCX)) ^c,oo)) 

PI^ + OO 

and such a travelling wave solution is called a soliton. We will be studying this type of 
solitary wave. In the heteroclinic case, the endstates are diherent, but must nonetheless 
satisfy a Rankine-Hugoniot-type condition. See [1] for more on these solutions, which are 
called kinks. 

From now on, we set c and Qc is a soliton, whose endstate, which we also set (thus 
we drop the index c), satishes 


Poodoi.Pco') ^ (^c« c) , (2) 

which means that Qoo is a saddle point for the hamiltonian ODE satished by Qc- Under 
this condition, we have that vanishes only once. 

The standard Lyapunov stability notion is that if solutions of equation ([T]) have initial 
conditions close to Qc, they remain close to Qdx — ct) at all times. But this notion is 
not satisfactory in describing the stability of travelling waves. Indeed, let c' 7 ^ c be close 
to c; we then have that the prohle Qc' is close to Qc, but, as the speeds are different, 
Qc'{x — c't) and Qc{x — ct) drift apart, despite their prohles remaining very similar. To 
see this, for a given t, compare Qc'{x — c't) with the translated prohle Qdx — c't) = 
Qdx — ct + {c — c')t). The correct notion of stability therefore stems from taking the 
difference of solutions with all the translated versions of Qdx — ct). A travelling wave 
solution will be considered stable if it is orbitally stable: for every e > 0 , there exists 
(5 > 0 such that if ||(po,wo) — Qc\\ < S, then 

sup inf \\{p{t),u{t)) — Qd' — a — cf)|| < e. 

t6K+ a6K 


The problem of orbital stability can be divided into two parts, depending on the type 
of perturbation we consider. ID perturbations are perturbations of Qc that depend only 
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on X and satisfy U 2 {t,x) = 0. The stability problem associated with these perturbations 
has been in part dealt with by S. Benzoni et al. in [1], and improved upon by Benzoni 
in [2] . A sufficient condition for orbital stability was obtained in the hrst paper using an 
argument by M. Grillakis, J. Shatah and W. Strauss [S], while the second article adds a 
sufficient condition for linear instability. See also J. Bowing mm for other stability 
results for the ID Euler-Korteweg system. 

The question of transverse stability deals with perturbations that also depend on the 
transverse variable y and have a 2D velocity held. So far, Benzoni in [2] and F. Rousset 
and N. Tzvetkov in [18] have proved linear instability. This occurs when the linearised 
equation around Qc has eigenvalues with positive real part. On one hand, Benzoni used 
Evans functions computations to get that orbitally stable solitons are transversally lin¬ 
early unstable. On the other hand, Rousset and Tzvetkov applied a simple, abstract 
criterion for instability in linearised PDEs with a hamiltonian structure in the case where 
the endstate of the soliton satishes (12|). We recall this criterion in Theorem 12.21 This 
criterion was applied to other equations with solitary waves in the same article: KP-I 
and Gross-Pitaevskii. 

The result of this paper is that the spectral instability mentioned above implies non¬ 
linear instability of Euler-Korteweg solitons. 

Theorem 1.1. Let Qc{x — ct) be a soliton solution to (QP such that the endstate Q^o 
satisfies Then there exist Sq and > 0 such that for every 0 < e < there exists 
an initial condition Uq = (po,uo) with 

\\Uo — Qc||j/s(]K2 ) < £ 

for some s > 0, such that, for every a G M, the solution U = (p, u) of (QP with this initial 
condition satisfies, at a time ~ ln(e“^), 

Moreover, the velocity u can be chosen to be potential: u = Vp. 

The proof relies on an argument originally by E. Grenier [7], in which one constructs 
an approximate solution G®'’’ to the equation based on a WKB expansion starting with 
the reference solution and unstable eigenmodes. An energy estimate on G — shows 
that, if there are enough terms in the expansion, this difference is small, thus, for times 
under the linear instability is dominant. Primarily used to obtain nonlinear insta¬ 
bility of boundary layers in numerous settings (unstable Euler shear flows and Prandtl 
layers [7|, Ekman layers for rotating fluids laiiaiia, Ekman-Hartman layers in MHD 
[5], Navier-Stokes with a boundary-layer-scale slip condition [I3]), the method relies on 
building approximate solutions from a compact set of unstable wave numbers, thus the 
idea has been transposed to showing transverse nonlinear instability of solitary waves, 
when these can be shown to be linearly unstable. F. Rousset and N. Tzvetkov have thus 
obtained nonlinear instability of solitary waves in many models: KP-I and NTS HZI, mul¬ 
tiple hamiltonian models including generalised KP-I and the Boussinesq equations [16] . 
and the free-surface water-waves equation [T9] . 
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Similarly to many situations, such as KP-I, our result extends to the |/-periodic frame¬ 
work M X Tl, where is the torus of length L > 0, in the following way: there exists 
a critical period Lo{c) > 0, depending on the soliton’s speed, such that if L > Lq{c), we 
have Theorem ll.il The proof is identical, and limited to above a critical period due to the 
loss of linear instability below this period (the set of unstable wavenumbers is bounded). 

Outline of the proof. The proof of Theorem 11.11 is in two parts. First, in section 
2, we build on Rousset and Tzvetkov’s linear instability theorem, deriving an important 
resolvent estimate for the Euler-Korteweg equation, linearised around Qc- This will allow 
us, in section 3, to build an approximate solution with the appropriate behaviour 
of being predominantly unstable for f ~ T^. An energy estimate on U — will then 
be used. This is particularly important ensure that the time of existence of U is large 
enough to get the desired amplihcation, as we only have local existence for solutions to 
o (see S. Benzoni, R. Danchin and S. Descombes in [3]). Article [3] also provides a blow¬ 
up criterion, but the instability phenomenon is not related to the blow-up if it occurs; 
indeed, the mechanism is also observed on systems that have global solutions (most of 
the previous examples, whether they concern boundary layer or solitary wave instability, 
fall in this category). Finally, combining the two will lead to the instability result. 


2 Linear analysis 

In this part, we consider the linearised problem about {pc,Uc) to obtain an essential 
ingredient in order to apply Grenier’s method. 

Considering that u is potential, we write the system on (p, ip), where u = 

f dtp + Vp ■ Vp + pAcp = 0 

\ dt^ + \\Vp>\^ = JF(p)Ap+lJG(p)|Vpp-po(p), 

which we linearise around {pc,Uc)- Having changed the space variable from x to x — ct 
(which turns the solitary wave into a stationary solution), we are interested in 

f dtp = {cdcc - Ucd^ - u'^)p - (p(.(9a; -h pcA)p . . 

\ dtip = {cd^-Ucd^)ip + {K{pc)A + K'{pc)p'^d:^-m)p, 

with m = K'{pc)p" + iF'(pc)(Pc)^ — 9oi.Pc)- We abbreviate the system by dehning two 
operators 

T 0 1 ^ 1 ^ ( -dx{K{p^)d:,) - K{p^)dly-m {uc-c)d^ \ 

V -1 0 y -5 x((mc-c)-) -dxipcdx) - pcdyy ) ' 

thus the system ([3]) can be written synthetically as dtV = JCV, where V = (p, v^). 


The hrst part of the linear analysis involves Ending unstable eigenmodes for ([3]). These 
are non-trivial solutions to the equation that can be written as V{t,x,y) = e'’"*e®^^n(x) 
for k ^ 0 and Re((T) > 0. We rewrite (jS]) using the Fourier transform on the transverse 
variable y: the equation dtV = JCV becomes av = JL{k)v with 


m 


-dx{K{pc)dx) + K{pc)k'^ - m {uc-c)dx A 
dx{^{uc c)’) dxi^Pcdx^ T pck' j 
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We begin by examining the existence of eigenmodes and the behaviour of a depending 
on k, and we follow up with an important resolvent estimate for JL{k). 

Proposition 2.1. Properties of the linearised equation. 


(a) The linearised equation is unstable, that is there exist eigenmodes written as 

V{t,x,y) = with v G and Re(cr) > 0, that solve For each k, 

the dimension of the subspace of unstable solutions of av = JL{k)v is at most 1. 

The instability is localised in the transverse Fourier space: there exists k^ax > 0 such 
that, for \k\ > fcmax, eigenvalues necessarily satisfy Re(cr) < 0. Let a{k) be the eigenvalue 
of JL{k) with highest real part. Then the function k i-G- Re((T(/c)) has a global maximum 
ao > 0 at a certain ko > 0. 

(b) IfV(t,x,y) = x), we define the following semi-norm for U: 

ll^(^)llx^ “ ll^l(^)ll_ffJ+l(R) + \\^xU2{t)\\^jQ^-^ + \k\ • 

It is essentially the norm of \k\U{t) plus the norm of U{t), omitting the 

norm of U 2 (t). 

Set 7 > (To, n > 0, s G M and let U solve 

dtU{t, x) = JL{k)U{t, x) + F{t, X, k), (4) 

with U (0, x) = 0, and, F satisfying, for every j < s, 

(5) 

uniformly for \k\ < /C ma x (Mg does not depend on k). Then U satisfies a similar estimate: 
for every j < s, and uniformly for \k\ < fcmax, we have 

\\dr^U{t)\\^,^<Cg{l+t)-^e'^h ( 6 ) 


Consequences of part (b): a quick energy estimate on the equation of U 2 yields 
that U 2 {t) G (as V9(0) = 0), and this If norm also satishes ([H])- We will therefore 

subsequently consider that the result is valid in if®, for any s > 0. 

By the Parseval equality, this result also implies identical if® bounds for (p, u) when 
written as (p, u) = (p, Vp) and (p, p) = f{y)e'^^^U{t, x) dk, with / G Indeed, 

norms of \k\‘^U can be replaced, using equation (jl]), by derivatives on x and t that satisfy 

dSD- 


2.1 Proof of Proposition 12.11 (a) 

Properties of eigenmodes 

2.1.1 Existence of unstable eigenmodes. 

The existence of unstable eigenmodes was shown by F. Rousset and N. Tzvetkov [18] 
using a general criterion for detecting transverse linear instability of solitary waves. We 
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have seen that equation ([3]) for functions written as V(t,x,y) = becomes an 

eigenvalue problem, that is 

av = JL{k)v (7) 

with J a skew-symmetric matrix and C a self-adjoint differential operator on 
whose domain is seen to be We have the following result for such systems. 

Theorem 2.2. (Rousset and Tzvetkov, /7^) 

If L has the following properties, 

• (HI) there exists fc mav > 0 and a > 0 such that L{k) > aid for \k\ > /Cmax/ 

• (H2) for every k 0, the essential spectrum of L{k) is included in [afe,-|-cxo[ with 

ak > 0 ; 

• (H3) L'[k) is a positive operator; 

• (Hf) the spectrum of L{0) consists of one isolated negative eigenvalue —A and a 
subset ofMH; 

then there exist a > 0 and k ^ 0 such that m has a non-trivial solution, and, for every 
unstable wavenumber k, such an eigenvalue a is unigue. 

This is shown by hnding k' > 0 such that L{k') has a one-dimensional kernel, and 
by using the Lyapunov-Schmidt method in the vicinity of this point; we do not detail 
the proof of this theorem. Proof that the linearised Euler-Korteweg system satishes the 
hypotheses of this theorem was also done in [TB], but we shall briefly recall this, as it 
contains some useful arguments for the subsequent points of Proposition 12.11 (a). 

(HI): using Young’s inequality, ab < |a^ -|- ^6^, with 5 = we quickly get that 


{L{k)v, v) > 


K{pc 


\dM"+{k^--]\v,\^ + Pc\dM"+le- 


K{pc) 


l^2p, (8) 


which is greater than a||17||^2 for \k\ large enough (remember that pc and K{pc) are 
positive). 

(H2): as \mix^±oo Qc{,x) = Qoo = (poo,Moo), with standard arguments [9], and using 
the fact that L{k) is self-adjoint, the essential spectrum of L{k) is given by that of 


.hoo (^) 


i^(Poo)(-<92^ + fc^)+^o(Po 
(o ^oo)dx 


[Un 


c)dx 


Poo(-5L + 


whose essential spectrum can be determined by using the Fourier transform in the x- 
variable and explicitly writing the eigenvalues p{f,k). We get that these are positive 
when k ^ 0. The essential spectrum of L{k) is equal to that of L^oik), so (H2) is verihed. 

(H3): we easily have L'{k) = dia.g{2kK{pc),2kpc). 

(H4): we apply the following lemma to h(0). 
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Lemma 2.3. Let L be a symmetric operator on a Hilbert space such that 


L = 



A 

L 2 


with L 2 invertible. Then, we can write 

{Lv, v) = ((L — AL 2 ^A*)vi, Wi) + (L 2 (t ’2 + L 2 ^A*vi), V 2 + L 2 ^A*vi^ . 


As a result, we write 

(L(0)u,u) = (Mui,ui) + f Pc 


dxV2 H - (mc - c)vi 

pc 


dx, 


with M = -dx{K{pc)dx-)-m- 


(Uc-C)^ 

Pc 


M is a second-order differential operator on which 


we can perform Sturm-Liouville analysis [B]. First, the essential spectrum of M is included 
in [a,-|-oo[ with a > 0; indeed M is a perturbation of M^o = ~K{poo)d1^ + doiPoo) — 
whose essential spectrum is positive under the assumption that Poo9oiPoo) > 
{uoo — c)^. Next, the function p^ is in the kernel of M, and it has one zero, so by Sturm- 
Liouville theory, M has a unique negative eigenvalue associated with an eigenfunction 
. Setting U 2 such that — c)n]", we have a generalised eigenfunction for 

L(0) (the second component is not in L^). By using approximations of U 2 , we see that 
(L(O)n, v) can be negative with U G H^, conhrming that L{0) has one negative eigenvalue. 


2.1.2 Localisation of instability and boundedness of unstable eigenvalues. 

We start be taking the real part of the scalar product of the eigenvalue equation ([^ 
by L{k)U\ we get that 'Re{a){L{k)U, U) = 0, as J is skew-symmetric. The operator L{k) 
satishes (HI) of Theorem 12.21 so, if Re((T) > 0, we must have 0 = {L{k)U, U) > a ||17||^2 
for \k\ > kma.x- Thus, the only function satisfying all = JL{k)U with Re((j) > 0 and \k\ 
large is f/ = 0; there are no unstable eigenfunctions for \k\ large. 

In order to get the boundedness of the unstable eigenvalues, we decompose L{k) as 
follows: L{k) = Lo(^) + -^1 with 

-dx{K{pc)dx) + K{pc)k'^ - mo 0 

0 -dx{pcdx) + Pck'^ 

where —mo(x) = max (—m(a;), |(7o(poo))- We compute the scalar product of ([7]) and 
LQ{k)v, and take the real part, which gives us 

Re{a){Lo{k)v,v) = Re{JLiv, LQ{k)v). (9) 

It is quickly noticed that there exists a > 0 such that 

Re{a){Lo{k)v,v) > aRe((T) (115,^^11^2 + 11^11^2 + \\vi\\l 2 ) ■ (10) 
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We shall now bound |Re(JLin, Lo(/c)n)| by the same norms as on the right. Note that 


f -dx{{uc-c)-) 0 \ 

^ V m-mo -{uc-c)dxj' 

Computing the scalar product directly, and bounding the terms involving and Uc in 
L°°, we get 

\Re{JLiv,Lo{k)v)\ < C [(1 + /c^)(|Re(ni, + |Re(9j,n2W2)|) + |Re(ni, 9j,n2)| 

+ 11*9x^111^2 + k"^ ||'^i||l 2 + 11*9x^211^2 
+ IK,* 9 xx'yi)| + \ivi,dl^V2)\ + \{\k\vi,\k\v2)\] 

The top line vanishes by integration by parts, and, also by integrating by parts, the hrst 
two terms in the hnal line are repeats of terms in the second line. Finally, by using 
Young’s inequality, the last term is bounded by ||n||^ 2 , thus proving that, combining 
with ([2]) and (ITU]) , there exists C > 0 such that 

Re(a) (119x^11^2 + k'^ 11^11^2 + 11^111^2) < C {\\d^v\\l2 + \\v\\l2 + \\vi\\l2) , 

which implies that Re((j) cannot be unbounded when positive. This ends the proof of 
part (a) of Proposition 12.11 


2.2 Proof of Proposition 12.11 (b) 

Resolvent estimate 

The proof of part (b) is split in two. First, we get the result for s = 0; we bound the 
norm of U by similar norms of F by using the Laplace transform and spectral arguments. 
The case s > 0 is then obtained by induction on s, the number of total derivatives (time 
and space). 


2.2.1 The case s = 0 


The proof of ([2]) for s = 0 relies on the Laplace transform, and is similar to the resolvent 
estimate proofs in irzi and na. Let (Jo < 7o < 7 . For f(t), we denote by /(r) the 
following Laplace transform; 


/(^):= 


'•+00 


exp(-(7o + ir)f)/(t) dt. 


Using the Laplace transform turns equation (j4]), dtU = JL{k)U + F, into an eigenvalue 
problem; 

(70 + ir)U{T) = JL{k)U{T) + F(r). (11) 

As 7 o > (Jo, 7o + iT is not in the spectrum of JL{k). Indeed, we can use the strategy 
employed to prove that hypothesis (H2) is satished to show that the essential spectrum 
of JL{k) is embedded in iM. Once again using the argument from [9j, we can examine 
the spectrum of the Fourier transform in x of JLoo{k), 


f -i{Uoo - C)i Poo(^^ + k"^) \ 

-iF(poo)(^^ + fc^) - fi'o(Pcx)) -i{uoo -c)i y ’ 
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which contains the solntions of the eqnation 

- 22^(Moo - c)X + pooK{p^){^‘^ + k‘^f + poo9o{poo)k'^ + {pooPoipoo) - (Uoo - = 0, 

which depend on {^,k). Using the positiveness of poo, K{poo) and condition ([2]), we get 
that the discriminant of this eqnation is negative for k) 7 ^ 0 , and clearly the only 
eigenvalne at k) = (0, 0) is zero, so the essential spectrnm of JL{k) is imaginary. 

As 7o + ir is not in the spectrnm of JL{k) for any r G M, the norm of the resolvent 
((70 + fT)Id — JL{k))~^ is nniformly bonnded for (r, k) in any compact snbset of It 
remains to show that, for \k\ < fcmax, there exists the following bonnd for |r| large. 

Lemma 2.4. If U solves /flil) . then there exists C, M > 0 such that, for |r| > M, 

||U(r)||^o<U||%)||^i. ( 12 ) 


Proof: we consider the scalar prodnct of the Laplace-transformed eqnation flTT]) with 
L{k)U, and write 

(70 + ir){L{k)U, U) = (F, L{k)U). (13) 

Note that {L{k)U,U) = {L{0)U,U) + K{pc)k'^\\Ui\\‘j ^2 -|-IU 2 I|^ 2 . Let ns concentrate 
on the term (L(O)U, U). Using Lemma [2.31 we know that it is eqnal to 


(L(O)U, U) = (MUi,Ui) + 


pc 


dxU 2 H- {uc - c)Ui 

Pc 


dx, 


The operator M, which we remind the reader is eqnal to —dx{K{pc)dx-) — m — ^ , 

and whose qnadratic form is dehned on has one simple negative eigenvalne, as well 
as a one-dimensional kernel containing p'^. 

Recall [/~ the generalised eigenfnnction corresponding to the negative eigenvalne of 
L(0). We do not have Uf G L'^, so we consider U~ = (f/f,0). We denote = (Pc,0), 
which is in the kernel of M 0 Id. We renormalise U~ and so that their norms are 
eqnal to 1. Let be orthogonal to U~ and we show that 


(L(0)U+,U+) >p\\U+ 


(14) 


for some 7 > 0. First, as Ui is not in the kernel or the negative eigenspace of M, we 
have {MUf,Uf) > a ihdih , since the essential spectrnm of M is inclnded in [a, -|-oo[. 
Thns, we already have 

(L(0)(/+,(/+)>a||(/+||7. (15) 

Bnt this does not snffice to get the Xq norm. Using the fact that pc is positive, there is 
a positive /3 snch that 


fp. 

Jr 

dxU+ + -{u,- c)U+ 
pc 

2 

dx > ft 

dxU+ + -{u,- c)Ut 
pc 





2C\{dxU+,Ut)\ 
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for some C > 0. We use Young’s inequality on the final term with 6 = (3j2C, thus there 
exists C" G M such that 

(L(0)f/+,f/+) > (a + C) ||(7+||^, + I . 

If perchance C is negative, we add x ffT5|) to the above, and obtain that there does 
indeed exist > 0 such that we have flT^ . and (L(/c)f/+, f/+) > r] 

We now write the orthogonal decomposition in of the hrst component, Ui = aU^ + 
hUi + Ui , and replace in ([13]). The eigenfunctions U~ and are hxed, so their 
norms are given constants. On the right-hand side, using integration by parts and basic 
estimates including Young’s inequality, we have 

|(F, mu)\ < I ||(7+||^. + CdIFII'l,, + a‘ + *2), 

k 

while on the left-hand side, we have 

(.mu,u) > ri(mu+\\l, + ||C/+||^j) - C(a^ + (|a| + |t|) ||f/+|ih). 

Taking the real part of flT^ and moving the negative part of the above to the right-hand 
side and once again applying Young’s inequality to absorb we get 

||(7+||i. <C(||f||y+o^+ (,=). (16) 

k 

To hnish off, we take the dot product of flTT]) with U~ and 17°. We quickly get 
(70 + *r)a =-(^,L(fc)J17-) + (F,17-) and (70 + *r)& =-([7, L(fc) J17°) + (F, t/°). 
Integrating by parts as usual, we get 

a^ + h^ < , ^1 |o (q' + ^' + 

7o + tI 

We see that if |r| is large enough, + 6 ^ can be absorbed by the left-hand side. Com¬ 
bining with flTHl) . we get the result. □ 


For the end of the proof of Proposition l2.ll (bh we start by using the Parseval equality 
in the following, 

/ e-"’»'||f/(()l|7 *= / Il7i)ll7*<c/ |7(()||„.<i( = C / \\F(t)\\„. dt 


Now recall the assumption on F, (Ej); we have 
f e-^^°^\\U{t)\\lo dt<CMo 


g2(7-7o)t 
) (1+^)“ 


dt < Co 


g2(7-7o)T 

(l + T)^' 


We inject this in the energy estimate on (|1|), that is 

|(||(7(0|li.)<C(||7(0|lt. + ||F(()|l?,.), 


and multiply the result by e integrate in time and we get the result. □ 
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2.2.2 The induction for s > 0 


The extension of Proposition 12.11 (b) to every s > 0 will follow the lines of the similar 
resnlt on the linearised water-waves eqnation in [T^]; it is done with a donble indnction, 
donble in the sense that one is embedded in the other. 

The hrst indnction is on s, the total nnmber of derivatives. Set s > 0, and we assnme 
that, for every s' < s and j < s', we have (Ej), that is 




K 


<a 


(1 


To get the wanted resnlt, we mnst prove that, for every 0 < j < s. 


2 

Xk 


:=\\drdiu\\i, + \k\^\\ardiu 


. 9 P 

r < c 

\l2 ^ 


2'yt 


(1 + 1 ) 


2n ' 


(17) 


where is the nsnal homogeneons Sobolev norm on M. This is done by indnction on the 
nnmber of space derivatives, j. Note that the norm (semi-norm if fc = 0) dehned here 
is a sort of homogeneons Sobolev norm expressed in the Fonrier space. Morally, at rank s 
of the indnction, we mnst get bonnds for the norms of terms involving s-|-l derivatives. 


Starting with j = 0, we are interested in the norm of dlU. Simply differentiate 
the eqnation, (jl]), s times with respect to t, and notice that W{t) = d^U{t) — d^U{0) 
satishes dfW = JL{k)W + G, with W\t=Q = 0 and ||G(t)||^s+i < 2Ms(l-|-we can 
re-nse the case s = 0 shown above. 


Now, let j > 0. To lighten the notations, we will write Vgj = d^~^dlV. Note that we 
want to control the norm of t/gj, which means s -|-1 derivatives in total, of which j + 1 
space or Fonrier derivatives. We apply d^~^di to the eqnation. This time, the derivatives 
do not commnte with JL{k), hence we consider 

dtUsj = JL{k)Us,j + J[di, L{k)]Us-j,o + Fsj := JMsjik)U + F^j. 


We take the real part of the scalar prodnct of this eqnation with Msj{k)U, which yields 
1 d 


2 dt 


{U,j, L{k)Usj) = -Re(f/,+i,„ [di, L{k)]Us-j,o) + Re(F,,„ Ms,,{k)U). 


To bonnd the second part of the right-hand side, we look more closely at the commntator 
term in M{k)U. We notice that there exist two sets of L°° matrices (m), m^)o<i<j+i snch 
that 

'i-i 


\,L{k)]dr^u = 


/i+i 

'^m]{x)Us-j+i,i ) + 
\i=0 




(18) 


, i=0 


We notice that all the terms, except the one with i = j + 1, have a total of s derivatives 
or less, and thns, nsing k"^ < fcmaxl^l, they are controlled by onr indnction hypothesis. 
Integrating the terms of {Fgj, L{k)Usj) involving j + 2 space derivatives and nsing as- 
snmption (E]) and Yonng’s ineqnality with a parameter rj to be chosen later, we obtain 
that the right-hand side is bonnded by 

fj 1-1 
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To deal with the hrst term of the right-hand side, we notice that Us+i,i has s -|- 1 deriva¬ 
tives, of which j space derivatives, hence the norm of falls under our induction 

hypothesis. We can thus use Young’s inequality and use decomposition flTSl) once again, 
and get 

1 d 

2'^^i^s,j',L{k)Us,j) < T] 

Finally, we integrate this in time, and recall ([ 8 ]) from the verihcation of the (HI) 
hypothesis of Theorem 12.21 which says that 

Us,) > i \\Us,\\\^ - c , 

in which the dual term can be moved to the right-hand side and controlled by the in¬ 
duction hypothesis. Note that we can take the norm of Ug-i^-i since the number of 
time derivatives is preserved (we do not encounter the second component of U which is 
not assumed to be in L^). In total, we therefore have 

\\U.j(T)\W < ||C/,,(t)|iy dt + Cj^^^. 

We choose r] in the Young inequalities above so that rj/r]' < 27 , and the Gronwall lemma 
gives us flT7|) for the couple (s, j). Both inductions are now complete. □ 

3 Nonlinear instability 

In this part, U = {p,u). Obtaining Theorem II . II relies on the construction of an approx¬ 
imate solution built around unstable eigenmodes of the linearised equation. In our 
case, this construction is classical and we will not write all the details of the calcula¬ 
tions (see also, for instance, IZIEI [T 6 | [T9]i. Energy estimates must then be obtained on 
U — to ensure that the approximate solution is close enough to the exact solution for 
long enough to see the difference between and Qc reach an amplitude 0{1). This is 
the more delicate part, as we remind the reader that only local existence is guaranteed 
for solutions of the Euler-Korteweg equation, thus we must also ensure that the solution 
U still exists when the instability appears. The closing argument is standard for this 
method. 

3.1 Construction and properties of the approximate solution 

For a whole number N independent of e to be chosen later, we will set 

N 

U""^ = Qc{t,x) + '^e^Uj{t,x,y), 
t=i 

with Uj potential: we set = (p'^P, Vj = {pj,ipj). This is expected to solve the 
Euler-Korteweg system leaving an error of order as follows, 

f dtp^^ + div (p^'PV^j’^P) = 

\ + i|V<p’^Pp + Po(p''P) = it:(p’^P)Ap’^P + |K:'(p’^P)|Vp’^Pp + 
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thus, expanding in these equations and isolating the terms of order e, we see that Vj 
solves the linearised Euler-Korteweg equation around Qc with a source term: 

dtVj = JCVj + Rj, (19) 

with Rj containing nonlinear interaction terms between the Vn with n < j, but with 
the sum on indices in each interaction term equal to j. For instance, we have -Ri = 0, 
and then, in the equation on ^ 2 , we see that -R 2 ,i = div (piV^Ji). In the equation on 
ifj, we notice that there will be quadratic terms similar to the one we have mentioned, 
but also cubic interaction terms, in the sense that they involve 14 ,i, 142 I 43 with 

ni + n 2 + Us = j. Likewise, R^'^ contains all the interaction terms whose sum of indices 
is greater than N. We do not detail the remainders any further. 


We now construct Vi as a wavepacket of unstable eigenmodes of the linearised equa¬ 
tion. Recall that fco > 0 is the global maximum for the function 

a : k ^ max{Re(A) | A G a{JL(k))}, 

where a{JL{k)) is the spectrum of the operator JL{k). We assume that this maximum is 
nondegenerate (although the method is easily adapted to the degenerate case, see [IS]). 
The function a is continuous, hence we can dehne 

Vi{t,x,y)= [ hik)R’^ye^^^K^ik,x) dk, 

Jr 


with fi{k) smooth, even, equal to 1 in the vicinity of k^ and supported in the set 
{k I d{k) > 3cro/4}, and w{t,x) = e'^^^Jvi{k,x) solving dtw = JL{k)w. For any s > 0, 
we consider the quantity 

dk = 

the second equality being the Parseval equality. We use the Laplace method around the 
critical point ko to get that ||Ri(t)||^s ~ so 

t—^+00 


s '=0 




s '=0 ' 


1 0 *^ 0 ^ 

Cl,,(1 + ty/* - + 


( 20 ) 


We get estimates on Vj by induction. Assume that / = supp(/i) is made up of two 
separate intervals around ±fco, and we set 

Vj{t, x,y) = J ''' J 4(^1) ■ ■ ■) x)e*^^^ • • • dki ■ ■ ■ dkj. 

Assuming that, for every n < j, 


\\vn{ki, - ■ ■ ,ki;t)\\jjs < Cnexp[n{d{ki) ^ - \-a{kn))t], (21) 

we get that Vj solves the linearised Fourier-transformed equation 

dtVj{ki, ■ ■■ ,kj) = JL{ki H-h kj)v{ki, ■■ ■ ,kj) + rj{ki, • • •, kj), 
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in which \\'rj{ki,---,kj)\\^^ < Cexp[n(a(/ci) + ■■■ + a{kj))t] by the structure of the 
remainder (combination of products between n„’s with n < j) and fl2T]) . Then, since, for 
k E I, a{k) > 3cro/4, the sum in the exponential is greater than cxo, and setting Vj\t=Q = 0, 
we can apply Proposition 12.11 (b) to get that Vj{ki, ■ ■ ■ ,kj;t) satisfies fl?T]) . We then use 
Parseval’s equality and the Taylor expansion of a around the critical point /cq to write 
that, for some (3 > 0, 


\\Vjikx,y)\\ls = 


ijkeji 


' fciH- \-kj=jk 


j{ki, • • •, kj] t, dki ■ ■ ■ dkj^i 


dk 


HO 


< 




' jk&jl 


'Em = l km=jk 




Integrate these gaussian functions (remembering that kj = jk — 
the desired inequality: for every j < 1, 

ll^'(^)llHqM2) (22) 

We now take a look at the remainder of the equation on which contains the 

interaction terms of the equation whose sum of indices is greater than N. Similarly to 
our proof of (12^ . we have 


3V 




j=N+l 


gjgicrot 

(1 +f)j/4’ 


(23) 


and, in what follows, we will be interested in times for which the smaller powers of 
e(l + are dominant. We set T* such that 


ee 


<^or* 


= K, 


(1 + T;)V4 

for 0 < K < 1 to be chosen later, and t = T* — r. Replace in fl23l) . and we have 

3N 


e^+^R^P(T; -r)| 


HO — 


max C, ) } 

,e{N+l,-,3N} J . 

N+1 —{N+l)aoT 






which, returning to the original time variable t, gives us, for t <T*, 

pN-\-l p{N-\-l)aot 

(24) 


3.2 Getting the instability 

If U is the solution of the Euler-Korteweg system ([T]) with the initial condition U (0) = 
17'^P(0), we will observe the instability by studying 

\\Uit) - QcimLHR^) ^ - ^'"WIIl2(R2) . 
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On one hand, we have — Qc = '^f=i and 


N 

i=i 


N 


> ||£f^l(^)||L2(]R2) 

1,2 (R2) i=2 


wu. 


J|IL2(R2) 


> r' 

- '(l+t)V4 


N 


i=2 


(1 + ty/^ 


by ([20]) and fl2^ . Taking times smaller than T*, we can consider that the snm on the 
right behaves like e^{l + and, replacing t by T* — r, we write 




r 


Il= > ''[c; 

> KC[e 




kC' 


_ '*'^2 „ — r7[jT\ 


Ci 


We notice that, for a given O > 0, there exists tc > 0 such that, for r > tc, 1 — Ce > 
1/2, so we set Ti > 0, independent of s, such that, for r > Ti, 


||(t/^p-a)(T; 


T 


2 ) > 




(25) 


On the other hand, we will use energy estimates to ensure that ||17(t) — l/'^^(t)||£, 2 ('R 2 ) 
is small. We will readily use those shown by S. Benzoni, R. Danchin and S. Descombes 
in [3], which are obtained by considering the equation on (G,z) = (G,u + iw), with G 
a primitive of the function p ^/K{p)/p and w = V(G(p)). The equation on 2 ; is a 
Schrodinger-type equation, written as 

dtz + u ■ Vz + iVz ■ w + iV(^(p)div z) = Q{p), 
with functions A and Q that we do not detail, while G satishes 

dtG + {u ■ V)G + ^(p)div u = 0. 


The approximate solution satisfies a similar system with a remainder term which we will 
denote 71 = {71g,7Zz)- From now on, we use tildes to designate the difference between 
the exact and approximate terms in this system, e.g. u = u — A = A{p) — A{p^'^). 
The difference (G, z) satishes the equation 

( dtG + u ■ VG + u ■ VG'^P + ^(p)div u + Rdiv = —TIq 
< dtZ + u ■ Vz + u ■ Vz^^ + iVz ■ w + iVz^'P ■ w (26) 

[ +iV(^(p)div z) + iV(^div z^p) = Q-n,. 

The energy estimates in on this equation involve multiplying by an adequately cho¬ 
sen gauge ^jJs■ In the potential case, this gauge is A{pY^‘^ (whereas in the non-potential 
case, a additional term is needed), and the weighted norm ||' 0 sA*-||^ 2 , in which is the 
standard Fourier multiplier for s derivatives, is equivalent to the standard norm. We 
send the reader to article [3] for details, in particular part 3 where the reasons for the 
simpler norm is explained, and part 6 where energy estimates on the difference between 


15 











a solution and a reference solution f fl26|) without the source term TZ) are obtained. 


Using the norms of Qc — Qoo as constants, the energy estimate reads 

^ (11^1 U2 + ||5||j^s)(1+||V2; •’||j^s+i)(H-|| V 2 ; ^11^^3-1 + 11511^^8)+^ ||7^||j^s. 

We rewrite this as follows: let fU = U — there exists a polynomial function Q such 
that, for s > 0 large enough (to handle the norms by Sobolev embedding), 

\\w{t)\\l. < rg(||u^p(^)-Qo.||^3 + ||u(^)||^3)||iu(^)||?,3 + ||i?^>’(^)||?,3 dz. 

Jo 

We will now choose N to get the right growth in time for W, as well as k to get the 
existence up to T* of the exact solution U. First of all, in the same way that we get fl25l) . 
we note that 


\\u^nt)-Qo.\\H. 


< \\Q-Qoo\\h- 

< \\Q-Qoo\\h- 


N 

+ 


i=i 

+ K, 


(1 + 


when t < T* — T 2 j with T 2 > Ti independent of e. We consider times t < Tw so that 
||hF(^)||j:/s < 1 and p(t,x,y) > 0 (no vacuum on the exact solution), and choose N so 
that, for t < Tw, 

2N(jq > g(||g — goo||ji/s +1 + k). 

A variant of the Gronwall inequality from na then provides us with 

(! + ,)(.,,)/. (27) 

for t < Tiy. Now, take t = T* — r: we notice that the right-hand side is smaller than 
which is therefore smaller than k if k < 1 is small enough. We now choose 
K, so that 2k < min pc, to ensure that there is no vacuum. We therefore have Tw > T* 
by a bootstrap argument. So, fl2711 is valid for t = T* — t with r > T 2 , and we have 

||W(T; - t)\\hs < (28) 


Recall that /i, which is the localising function for the unstable term f/i, is supported 
in I which is made up of two closed intervals that do not contain 0. Set / a smooth 
function such that f{k) = 1 for k E I and / = 0 in the vicinity of 0, and let us dehne fl, 
a Fourier projector on frequencies in /, by 


I'y(Jlu){x,k) = f{k){J^yu){x, k). 

As f/(0,a:, p) = Qc{x)-\-eUi{Q,x,y), nf/|t=o = Uq. Moreover, for any a G M, the difference 
Qc{x — ct — a)— Qc{x — ct) does not depend on y, hence n(gc(- — ct — a) — Qc(- — ct)) = 0. 
We can now combine fl25H and fl28|) to show the instability: we have, for any a. 


\\U(t) — Qc{- — a — ct)\\j^2^T^2) \t=T*-T > 

> 


||n(u(f)-Q,(--cf)|| L2(R2) U=T*-t 


kC[ 


e-^or 


1 


2k^C'o 


^-Naor 
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For r > Ts > r2, we have the last exponential on the right smaller than 1/2, and, as a 
result, letting t' > T3 be fixed, independent of e, 

kC 

:=5o. 

The number we have found is positive and does not depend on e: Theorem 11.11 is 
proved. □ 
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